Tsallis statistics generalization of non-equilibrium work relations 
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Introduction- Recent advances in theory of non- 
equilibrium statistical mechanics have established the 
methods to calculate the free energy differences between 
the two equilibrium states of the driven system from the 
non-equilibrium work measurements [lMSI]- These meth- 
ods are generally called as non-equilibrium work rela- 
tions and in particular named as the Jarzynski equality 
[H and the Crooks work fluctuation theorem Q. Con- 
sider a system initially in equilibrium at temperature (in- 
verse) (3 — 1/kT (k is the Boltzmann constant), which 
is externally driven from its initial equilibrium state A 
to final equilibrium state B by non-equilibrium process. 
Let Pf[j f ] be the probability of the phase space trajec- 
tory 7 F , for the system driven between the two states in 
forward direction. This satisfies the Crooks work fluctu- 
ation theorem @, 0] , 



P F [l F ] 



a l3(W-AF) 



(1) 



where W is the work performed on the driven system, 
AF is the free energy difference between the two equi- 
librium states and Pr[j r ] is the probability of the phase 
space trajectory j , for the system driven in the reversed 
direction. This is the direct relation between the work 
dissipation and the ratio of probabilities for the forward 
and the reversed trajectories and its integrated version is 
Jarzynski equality 



aiccto 
0,1, 



(exp[-/3W}) =exp[-£AF]. 



(2) 



The average (...) is over a statistical ensemble of realiza- 
tions of a given thermodynamic process. 

In Crooks work fluctuation theorem, the probabili- 
ties of non-equilibrium forward and reversed trajecto- 
ries are related by taking the initial conditions from the 
Boltzmann-Gibb's (BG) equilibrium distribution. There 
are very few studies on work fluctuation theorem relat- 
ing the non-equilibrium forward and reversed trajectories 
taken from other statistical distribution [a, 0] ■ I n partic- 
ular, finite bath work fluctuation theorem has been de- 
rived earlier [3] , which includes the microcanonical work 
fluctuation theorem and the Crooks work fluctuation the- 
orem as the two limiting cases. This has been obtained 



by considering Tsallis statistics [8] as a finite heat bath 
statistics Q- However, this theorem generally contains 
two temperatures instead of one as observed usually in 
non-equilibrium work relations (Eqs[T]and[2]). There has 
not been any generalized connection established between 
Tsallis statistics and non-equilibrium work relations at 
a single temperature. Since the generalized connection 
should exhibit interesting applications in complex sys- 
tems 0, in this paper, we derive the g-statistics gen- 
eralization of Jarzynski equality and the Crooks work 
fluctuation theorem for the classical system driven be- 
tween two equilibrium states by a non-equilibrium pro- 
cess using Tsallis statistics. 

The theory of Tsallis statistics based on a generalized 
form of entropy, S q , characterized by the index g 6 i?, 
such that q — 1 recovers the standard theory of Boltz- 
mann and Gibbs. This generalized (Tsallis) entropy is 
given by the expression [S] 



1 



(9-1) 



(3) 



where k is the positive (Boltzmann) constant and 

w 

c = £# (4) 

i=l 

Here, w is the total number of microstates of the system 
and Pi is the probability of the system at microstate i. 
In the limit q — > 1 one can recover BG entropy 



'BG 



k^Pi In Pi 



(5) 



Preserving the standard variational principle, Tsallis es- 
tablished the canonical generalized distributions and its 
refinements @, [llj] . Using Tsallis statistics, Chame and 
Mello have derived the generalization of the fluctuation 
dissipation theorem 12|. Tsallis nonextensive statisti- 



cal mechanics is also considered as an approach to non- 
equilibrium stationary states of small or complex systems 
0]. However, its equilibrium formulation remains to be 
valid for obtaining thermodynamic properties of the equi- 
librium system (jj EH • There are four versions of Tsallis 



2 



statistics 14|, in particular, we use the most widely ac- 
cepted third constraint (escorted probability) [§] formu- 
lation of Tsallis statistics. 

In Tsallis third constraint formulation, the generalized 
equilibrium canonical distribution at ft is given by 11 1 



!-(!-?) 



fth - u g ] 



z n 



expj-/3(ej - U q )/c q ] 



(G) 



where e; is the energy of the ith microstate, U q is the 
normalized constrained of internal energy which is given 
by [HI, [l5| 



1 10 



(7) 



and Z q is the ^-generalized partition function which is 
given by Q 



i-(i-g) 



fth - u q 



(8) 



The normalization condition of pi leads to the relation 



[HI 



-q — ^ q 



(9) 



This modified formalism also becomes ordinary canonical 
ensemble theory in the limit q — > 1 [9| with c g= i = 1 and 

w 

Z q=1 = exp[ftU] J^exp [-fo] = exp [/3(f/ - F)], (10) 



^ fiiexp [-/3ei] 



where the internal energy 

u = expjftU] 
Z q =i 

e t exp [-fkj\ 

XT eX P \-P e i] 

and the free energy 

w 

F = -kT in exp [-fte,] . 



(11) 



(12) 



Consider a system in an initial macrostate A (for exam- 
ple closed system of volume Vi) which is in equilibrium 
at ft. The probability for the system in a microscopic 
phase-space (microstate) Ta is given by 111. Il6| 



p(r. 



Z q (A) 



ft[H{T A ) - U q {A)] 
c q (A) 



(A3) 



where H(Ta) is the Hamiltonian for the system in a mi- 
crostate T A , U q (A) is the internal energy which is the 



(escorted probability) weighted Hamiltonian eigenvalue 
11] averaged over all microstate in an initial equilibrium 
state A (see, Eqs. HI [7] and [5]) and 



Z q (A) = 



l-(l-«) 



ft[H(T A ) - U q (A)} 



c q {A) 



with 



(14) 



(15) 



Suppose the given system evolves in time under Hamil- 
tonian dynamics and reaches a different macrostate B 
which is to be in equilibrium at same ft. The probability 
distribution for the system in a microstate is given 
by, 



1 



Z q (B) 



1 - (1 - q) 



ft[H(T B )-U q (B)} 
c q {B) 



where H(Tb) is the Hamiltonian for the system in a mi- 
crostate r B , U q (B), Z q (B) and c q (B) have same meaning 
as above but for the macrostate B. In order to derive the 
non-equilibrium work relations for a given ft, we formu- 
late the problem as follows. 

Setup- Consider a classical Hamiltonian system in a 
macrostate A which is initially in equilibrium with reser- 
voir at inverse temperature ft. Let At be an external 
protocol applied in the arbitrary time interval r to drive 
the system from its initial equilibrium state A to an- 
other state B at constant bath temperature ft. It is as- 
sumed that the final state B is not necessarily to be in 
equilibrium. However, the system in the state B at the 
constant bath temperature relax towards the equilibrium 
state B for the same ft without doing any work [5j. Let 
H(T t ,Xt) is the Hamiltonian with externally controlled 
time-dependent protocol At and the phase-space coordi- 
nates of the system, L t at a particular time t. At t = 0, 
the system Hamiltonian which is in any one of the mi- 
crostate is H(T ,X ) = H(Ta); and at time t = r, 
the system Hamiltonian is H(T T ,A T ). Let 7 denote the 
entire trajectory of the driven system from t = to r. 
One can obtain the statistical ensemble of possible real- 
izations by performing the above process repeatedly. In 
following the refs [H, la 0] , the work performed on the 
system for a given trajectory can be defined as 



W = H{T T , A T ) — H(Tq, \q) 

= H{r T ,x T )-H{r A ). 



(17) 



It should be noted that the work defined in Eq. (fT7]l is 
different from the o-dependent work as given in ref.[llj 
(see also refs. 3, 17jj). Different definition of work and its 
physical meaning has been discussed in detail in refs. [4, 

MM. 
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During the time interval r in which the system is 
driven, we assume that the reservoir should always be 
in equilibrium at a given /? [2(|. In such a case, the total 
heat transferred by the system can be written into two 
part as 



Q = Qt + Q r , 



(18) 



where Q q is the g-dependent heat transferred between 
the system and the reservoir which should preserve the 
(g-dependent) equilibrium nature of the reservoir and Q r 
is the heat transferred when the system relaxes towards 
the required equilibrium state from the final state B at a 
given j3 2l|. There is no work performed on the system 
during relaxation We can define the heat transfer 
due to relaxation as 



H(T B )-H(T T ,Xr). 



(19) 



Since the energy conservation 
equilibrium process [I 



22, 2 



is also valid for non- 
, the above driven system 
should obey the principle of energy conservation for any 
microscopic trajectory as 



W + Q = AU, 



(20) 



where AU q — U q (B) — U q (A). The principle of energy 
conservation as given in Eq. (|20p for the non-equilibrium 
process [23| is different from the g-generalization of first 
law as proposed in refs.fll, Til 17 [ for equilibrium system. 
Therefore, the g-dependent heat transferred by the driven 
system between the two equilibrium state is 



Qq = AU q -W-Q r 



(21) 



= AU q — [H(Tb) — H(Ta)} 



which depends only on the initial and final system states. 
Since Q q is independent of the non-equilibrium trajecto- 
ries of the driven system, we impose the condition for 
the system relaxing towards the required (g-dependent) 
equilibrium state at a given j3 as 



Cq(B)=C q (A) 



!-(!-?) 



PQ d q 

c q (A) 



(22) 



Using Eq.Q, the above equation can be written as 

l-Cg(g) 1 - Cg(A) 



PQ a a 



(23) 



-ISq(B) - S q (A)] 



The above equation provides the consistent usage of j3 = 
1/kT as the (inverse) reservoir temperature in Tsallis 
statistics [l7j |. 

The usage of thermodynamic laws are not mandatory 
for the proof of non-equilibrium work relations 0, Hj]. 
Since both internal energy, applied work and heat are for- 
mulated clearly for the above driven system [2^|, we use 
Eq. ([20|) and Eq. (f22|) and obtain the ^-generalized non- 
equilibrium work relations in which the initial probability 
distributions (Eqs. [TBI and fTBj) are taken from the Tsallis 
statistics. 

q- generalized Jarzynski equality- In order to obtain the 
g-generalized version of Jarzynski equality for the above 
driven process, one can take the following g-exponential 
average over an ensemble of realizations in which the ini- 
tial distribution is taken from Tsallis statistics. 



-l3(W + Q T -AU q ) 



-l3(W + Q r -&U q ) 



P(T A ) dT A .{25) 



Using the g-exponential identity e x q e v q = eq x+ l' + ( 1 ~ q ' xy ' 
9} and Eq. (fT3| . the integral of the above equation can 
be rewritten as 



-l3(W + Q r -&U q ) 



where, 



/3[W + Q r - AUq] 
c g (B) 



(27) 



and 



JV„ = - 



g^ [l + d-,)MJ. (28) 



(29) 



Using Eq.(|20|) and Eq. (|22|) . one can rewrite 
1 + -q)M q = l + (l-g) ^ "• 



c q (B) 
= 1 | c q (A)-c q (B) 



c q (B) 



c q (B)' 



Using Eqs. ([T7l and [T9l) . one can obtained 



where a q = S q (B) — S q (A) is the g-dependent change in 
entropy of the equilibrium system From the above 
condition (Eq. |2"2"|) . the g-dependent change in (equilib- 
rium) reservoir entropy for the driven non-equilibrium 
process is obtained as 23h25j 



as: 



(24) 



Mq + N q = -JL[H(T B )-U q (B)]. (30) 

Cq{B) 



Therefore, Eq.pfJ]) becomes 



-P(W+Q r -AU q ) 

, CqTB) 



Z q {A) 



0lH(r B )-U q (B)] 
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Since Hamiltonian dynamics preserve the phase-space 
volume, cLTa = dTs 0,01 an d using Eq. fTB]) we can 
rewrite the above equation as 



-P(W + Q r -AU q ) 
<=<i( B ) 



/P ( r 



B )dT B (31) 



Z g (B) 
Z q (B) 

z q (Ay 



We have obtained q-generalized version of one of the non- 
equilibrium work relation. It should be noted that /3 ap- 
peared in the g-exponential average to ensure that the 
temperature of the reservoir remains same, however, one 
does not know anything about system temperature dur- 
ing the driven process. Further, Q r and c q {B) in the 
above relation also takes care of the heat exchange due 
to relaxation of the system towards the final equilibrium 
state [2l| . This may provide the possible physical mean- 
ing of the above average for the driven non-equilibrium 
process instead of thinking as an adhoc method 0, [2l| . 

In the limit q — >■ 1, exp g (x) = exp(x) and c q = 1 9], 
then using Eqs. (|10H12|) . Eq.fl3T]) becomes 

(cxp[-f3(W + Q r - AU)}) = exp [(3(AU ~ AF){32) 
(exp[-(3(W + Q r )}) = exp[-/3AF], 

where AU — U(B) — U(A) is the change in internal en- 
ergy and AF = F{B) — F(A) is the change in equi- 
librium free energy of the BG canonical system. Thus, 
we can obtain more general form of Jarzynski equality 
which includes heat due to system relaxation [2l| for the 
BG canonical system in the limit q — > 1. If Q r is within 
the measurements/numerical error for work calculation 
in experiments/simulations, the original Jarzynski equal- 
ity (without the heat term due to relaxation 0, Hlj) can 
provide the reliable estimates of the free energy differ- 
ences. 

q- generalized Crooks Work fluctuation theorem- In or- 
der to derive the q-generalized version of the Crooks work 
fluctuation theorem, we proceed with the problem anal- 
ogous to ref.[H as follows. Since H(T , Ao) = H(Ta), the 
probability of the phase space trajectory -f F , for the sys- 
tem driven between the two equilibrium states obtained 
from the initial equilibrium distribution in the forward 
direction (A to B) is given as Q 



Pph 1 



(33) 



Z g (A) 



l-(l-9) 



P[H(T A ) - U q (A)} 
c q (A) 



Suppose the system is driven from equilibrium state B 
to state A using the time reversed protocol, Af = Af_ t , 
H(T R ,X R ) is the Hamiltonian for the externally con- 
trolled time-dependent protocol A^ and, T R is the phase- 
space coordinate of the system at a particular time t. At 
t = 0, the system Hamiltonian which is in any one of the 



microstate Tb is H(T R , X R ') = H(Tb)' : and at time t = r, 
the system Hamiltonian is H(T R , \ R ). The probability 
of the phase space trajectory 7^, for the system driven 
in reverse direction obtained from the initial equilibrium 
distribution is given as Q 



Pr[i r ] = P e B q {K) 



(34) 



Z q {B) 



l-(l-ff) 



P[H(T B ) - U q {B)\ 



c q (B) 



Using Eq.(H2]), Eq $19$ and Eq.®, Eq.(|33]) can be 
rewritten as 



Pph 1 



Ppb 1 



1 



Z q {A) 



1 



i-(i-s) 



P[H{T B )-U q {B)+> 



c q {A) 



Z q (A) 
1-(1- ,) 



l-(l-9) 



PQ 



Using Eq. 



[ J Z q (A) 



c q (A)_ 

P[H(T B ) - U q (B)} 
c q (A)[l-(l-q)^ 

and Eq.flM]), Eq.flSS]) becomes, 

PQt 



(35) 



l-(l-«) 



c q (A) 



Pr[i H }- (36) 



We can get from Eq.(gD]) that 



P F [l F ] 
Pr\i r ] 



Z g (B) 
Z g (A) 
Z q (B) 
Z q (A) 



l + (l-«) 
P(W 



p [W + Q r - AU q ] 



exp„ 



c q (A) 
Q r - AU q ) 



l-q 

(37) 



c q (A) 



We have obtained the g-generalized version of another 
non-equilibrium work relation. In the limit q — > 1, 



exp q (x) — exp(x) and 
[T2J) , Eq.® becomes 

Pf[ 7 Fi 



then using Eqs.dTDl- 



Pr[j r ] 



= exp [p(W + Q r - AF)} . 



(38) 



Thus, we can obtain more general form of Crooks work 
fluctuation relation which includes heat due to system 
relaxation 21] for the BG canonical system in the limit 
q^l. 

In order to obtain the g-generalized version of Jarzyn- 
ski equality from the g-generalized work fluctuation re- 
lation, one can take the following q-exponential average 
over ensemble of realization in forward direction as 



-/3(W + Q r -AU q ) 



-P(W + Q r -AU q ) 
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Using Eq. ([37|) . the above equation can be rewritten as 



-t3(W + Q r -&U q ) ' 



Z g (B) T 
Z q (A) «' 



(40) 



where 



-P(W + Q r -AU q ) fl(W + Q r -AU q ) 

I q — I V q tq 



P R [l R ]d lR . (41) 



Since Hamiltonian dynamics preserve the phase-space 
volume, d'jF = dj R 0, Using the ^-exponential iden- 
tity e x q e v q = e]? +y+ ( 1 ~ q ' xy ] Q : W e rewrite the integral of 
the above equation as 



where 



D q = 



1 



J(W+Q r 



- AU ^Pnh R ] d lR . 



1 



c q (A) c q (B) 
Using Eq.flU and Eq. 

1 1 



-(1-9) 



P(W + Q r - AU q 



(42) 



(.43) 



Dq = 



Cq(A) Cq (B) 

the above equation becomes, 
[c q (B)-c q (A)Y 



[Cq(A) Cq(B) Cq{A)Cq{B) 



0. (44) 



Since e° = 1, I q = J Pr[j r ] d lR = 1 and Eq.® be- 
comes, 



exp g 



-p(W + Q r -AU q ) 



c q (B) 



Z q {B) 

z q {Ay 



We have obtained the ^-generalized version of Jarzynski 
equality. 

Conclusion- We have derived the more general form of 
Jarzynski equality and the Crooks work fluctuation the- 
orem which includes the heat due to system relaxation 
in the framework of Tsallis statistics. Our gen eral result 
may resolve the criticism raised earlier (EL l2lj for origi- 
nal Jarzynski equality. In Tsallis third constraint formu- 
lation, one cannot directly obtain the canonical proba- 
bility distribution because the distribution (Eq|6]) is self 
referential 0,[I3|- Since pi depends upon c q , one should 
iterate Eqs.([51 15] and O repeatedly until numerical con- 
sistency is achieved. Since Tsallis distribution provides 
potential application in various complex systems, we have 
utilized the self referential nature of the Tsallis distribu- 
tion and have obtained the g-generalized version of non- 
equilibrium work relations. 

There is a general impression among few of us that 
Tsallis's formalism has nothing to do with equilibrium 
statistical mechanics. If the system relaxes towards the 
equilibrium, one may not rule out the equilibrium formu- 
lation of the Tsallis statistics. Although we have taken 
the initial distribution as the equilibrium, our formula- 
tion may also be applicable for non-equilibrium station- 
ary state conditions. 
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